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SOME COMBINATORIAL ASPECTS OF QUANTUM FIELD 

THEORY 

ADRIAN TANASA 

Abstract. In this survey we present the appearance of some combinatorial notions 
in quantum field theory. We first focus on graph polynomials (the Tutte polynomial 
and its multivariate version) and their relation with the parametric representation of 
the commutative $ field theory. We then generalize this to ribbon graphs and present 
the relation of the Bollobas-Riordan polynomial with the parametric representation 
of some $ 4 field theory on the non-commutative Moyal space. We also review the 
role played by the Connes-Kreimer Hopf algebra as the combinatorial backbone of 
the renormalization process in field theories. We then show how this generalizes to 
the scalar $ 4 field theory implemented on the non-commutative Moyal space. Finally, 
some perspectives for the further generalization of these tools to quantum gravity 
tensor models are briefly sketched. 



1. Introduction 

Within the framework of quantum field theoretical models, it has been shown recently 
that the role of various combinatorial notions is a crucial one. 

Inside the area of theoretical physics, quantum field theory (QFT) represents an im- 
portant part; its mathematical formalism was proven successful not only in elementary 
particle physics (the celebrated Standard Model) but also in condensed matter physics. 

Like in other domains of theoretical physics (exact solutions of statistical-mechanical 
problems, random matrix models, integrable systems and so on), also in QFT a key 
role is played by combinatorics. Thus, perturbative QFT relies on graph theory, namely 
on (Feynman) graphs which appear in the expansion with their combinatorial weights; 
moreover, as we will present in this survey, graph polynomials can be naturally related 
to polynomials appearing in some representation of the Feynman amplitudes in QFT. 
Furthermore, renormalization, which lies at the heart of QFT, can be described through 
some appropriate combinatorial Hopf algebras, as we will see in the sequel. 

On the other hand, techniques of analytic combinatorics (such as the Mellin trans- 
form or the saddle point method) are also used in QFT for computing and analyzing 
quantities like Feynman integrals (which are associated to Feynman graphs). 

It is worth emphasizing that, in the last years, both combinatorialists and theoretical 
physicist in general (or field theorists in particular) became more and more aware of this 
strong relation between their domains. This is a natural tendency, since the unfolding of 
new ideas in physics is often tied with the development of new combinatorial methods. 
We can thus speak nowadays of an emerging domain of combinatorial physics. 



As already announced above, in this survey we focus on the appearance of combi- 
natorial notions in QFT, such as graph and map polynomials and combinatorial Hopf 
algebras. 
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Thus, in graph theory the celebrated Tutte polynomial (see (72J E3]) is known to 
characterize in a particularly elegant way a generic graph. On the other hand, the 
Feynman integral of such a graph can be represented, through the parametric represen- 
tations, using some polynomials in a set of parameters associated to the edges of the 
graph. In this survey we present the relation between the Tutte polynomial (or, more 
exactly, its multivariate formulation (65]) and the polynomials of this QFT parametric 
representation [Tj [TBI |4~4"]. 

On the other hand, we also present how the Connes-Kreimer Hopf algebra of Feynman 
graphs [21] encodes in a powerful manner the combinatorial structure of perturbative 
renormalization in QFT. 

We then generalize this framework and lift it to ribbon graphs (also known as maps, 
see for example [18] and references therein). At this level, the Tutte polynomial can be 
replaced in a natural manner by the Bollobas-Riordan polynomial [9j [10]. From the 
QFT point of view, the appropriate models are the non- commutative ones (we present 
here the scalar $* 4 model implemented on the Moyal space). In this case also, a 
parametric representation generalizing in a highly non-trivial manner the commutative 
one can be defined [HJ [68] . One can then relate the Bollobas-Riordan polynomial to 
the polynomials of the parametric representations of these QFT models on the non- 
commutative Moyal space , r FI] . 

Furthermore, a ribbon graph version of the Connes-Kreimer Hopf algebra of Feynman 
graphs can be defined [TO"! ITT] . This algebra is proven to lie behind the renormalization 
of QFT models on the non-commutative Moyal space. 

2. Some algebra 

In this section we briefly recall the definitions of the algebraic notions that will be 
used in the sequel. 

Definition 2.1 (Algebra). A unital associative algebra A over a field IK is a K-linear 
space endowed with two algebra homomorphisms: 

• a product m : A <S> A — > A satisfying the associativity condition 

mo (m® id)(r) = mo (id®m)(r), for all T G A® 3 ; (2.1) 

• a unit u : K — > A satisfying 

mo {u® id)(l <g> T) = T = mo (id® u)(V ® 1), for all V 6 A. (2.2) 

Definition 2.2. A (coassociative, counital) coalgebra C over a field IK is a K-linear 
space endowed with two linear homomorphisms: 

• a coproduct A : C — > C Cg> C satisfying the coassociativity condition 

(A® id) o A(r) = (id® A) o A(r), foralireC; (2.3) 

• a counit e : C — > K satisfying 

(e®id) o A(r) = T = (id®e) o A(r), for all T e C. (2.4) 

Definition 2.3. A bialgebra B over a field K is a K-linear space endowed with both an 
algebra and a coalgebra structure (see Definitions 12.11 and 12.21) such that the coproduct 
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and the counit are unital algebra homomorphisms (or, equivalently, the product and 
unit are coalgebra homomorphisms): 

A o mg = m mB o (A <g> A), A(l s ) = \ B ® 1b, (2.5a) 

£ o m B = m K o (e ® e), e(l B ) = 1. (2.5b) 

Definition 2.4. A graded bialgebra is a bialgebra graded as a linear space, 

oo 

£ = 0£ (n) , (2.6) 

71=0 

such that the grading is compatible with the algebra and coalgebra structures: 

n 

B (n) B (m) £ g(n+m) and A g(n) g Q g(fc) g, g(n-A). (2.7) 

fc=0 

Definition 2.5. A connected bialgebra is a graded bialgebra jB for which i^ ) = w(K). 

Definition 2.6. A Hopf algebra % over a field K is a bialgebra over IK equipped with 
an antipode map S : H —> H obeying 

m o (S <g) id) o A = u o e = m o (id ® S 1 ) o A. (2.8) 

A practical introduction to Hopf algebras for a combinatorial physicist can be found 
in [26]. For further details on this topic, the interested reader is referred to [41 J or [24J, 
for example. 

3. Graph theory — the Tutte polynomial 

3.1. Some notions of graph theory. In graph theory, one has the following definition 
(see, for example, [74]). 

Definition 3.1. A pseudo-graph T is defined as a set of vertices V and a set of edges 
E together with an incidence relation between them. 

This means that multiple edges (i.e., edges connecting the same two vertices) and 
(self-)loops (i.e., edges connecting a vertex to itself) are allowed. Nevertheless, we will 
simply refer to these pseudo-graphs as graphs in the rest of this survey. 

One can extend the definition above such that a distinct type of edge — the external 
edge — is permitted. Such an edge is attached to only one vertex. 

In general, the terminologies used by graph theorists or by field theorists are different. 
For the sake of completeness, we present in this section both of them. Nevertheless, in 
the rest of this survey we use the graph theorists' language. 

Definition 3.2. 

(1) The number of edges at a vertex is called the degree of the respective vertex 
(field theorists refer to this as the coordination number of the respective vertex). 

(2) An edge whose removal increases the number of connected components of the 
respective graph is called a bridge (field theorists refer to this as a 1-particle 
reducible edge). 

(3) A connected subset of edges and vertices, the number of edges being the same 
as the number of vertices, which cannot be disconnected by removing any of the 
edges is called a cycle (field theorists refer to this as a loop). 
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(4) A (self-)loop (see above) is called a tadpole edge in QFT terminology. 

(5) An edge which is neither a bridge nor a self-loop is called regular. 

(6) An edge which is not a self-loop is called semi-regular. 

(7) A graph with no cycles is called a forest. 

(8) A connected forest is called a tree. A spanning tree is a tree connecting all the 
vertices of the graph. 

(9) A spanning forest with two connected components is called (in QFT terminol- 
ogy) a two-tree. 

(10) The rank of a subgraph A is defined by 

r(A) := \V\ - k{A), (3.1) 

where k(A) is the number of connected components of the subgraph A. 

(11) The nullity (or cyclomatic number) of a subgraph A is defined by 

n{A) :=\A\-r{A). (3.2) 
We illustrate this notion by the example of Figure [TJ 




Figure 1. An example of a graph (with seven edges and six external 
edges). We chose a spanning tree and we label its edges by 1,2,3,4. 
The set {1,2,4} is a two-tree; one has two connected components (the 
first one formed by the edges 1 and 2 and the second one formed by the 
edge 4). The external edges are attached to one of these two connected 
components. 

3.2. Graph polynomials and the (multivariate) Tutte polynomial. One can 

define two natural operations for an arbitrary edge e of some graph T: 

• the deletion, which leads to a graph denoted by T — e, 

• the contraction, which leads to a graph denoted by T/e. This operation identifies 
the two vertices V\ and at the ends of e into a new vertex t> 12 , attributing all 
the edges attached to V\ and Vi to V\2, and then it removes e. 

Remark 3.3. If e is a self-loop, then T/e is the same as T — e. 

For an illustration of these two operations, see the example in Figure [2j 

Let us now give a first definition of the Tutte polynomial. 
Definition 3.4. If T is a graph, then its Tutte polynomial Tr(x,y) is defined by 

T T {x,y) := Y,{x-l) r{E) - r{A \y-V n{A) - (3-3) 

AcE 
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Figure 2. The deletion/contraction of some graph. One is left with 
various possibilities (here five) of terminal forms (that is, graphs with 
only bridges or self-loops). 



A fundamental property of the Tutte polynomial is the deletion/contraction relation. 

Theorem 3.5. If T is a graph, and e is a regular edge, then 

Tr(x, y) = T T/e (x, y) + T r „ e (x, y). (3.4) 

Let us remark that this property of the Tutte polynomial is often used as its definition. 
However, for this property to become a mathematical definition, one needs to complete 
it by giving the form of the Tutte polynomial on terminal forms. Namely, for graphs 
with only irregular edges, m bridges and n self-loops, the Tutte polynomial is given by 

T T {x,y):=x m y n . (3.5) 

A multivariate version of the Tutte polynomial exists in the literature [SS] (see also 
[28, 29J). The main idea is that, instead of having a single variable, y. for the number 
of edges, one introduces a set of variables /3 2 , • • • , P\e\, one for each edge. This leads 
to the following definition. 

Definition 3.6. If T is a graph, then its multivariate Tutte polynomial is defined by 

Z r (q,{(3}):=J2^ A) U^- ( 3 - 6 ) 

AcE e£A 

Similarly, one can prove that the multivariate Tutte polynomial (13.61) satisfies a dele- 
tion/contraction relation, for any edge e. The definition of the polynomial on terminal 
forms (graphs with v isolated vertices) is: 

Z r (q,{(3}):=q v . (3.7) 
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Through direct inspection, one can prove the following relation between the Tutte 
polynomial ( 13 .3p and its multivariate counterpart (13. 6ft : 

[q- v Z T {q^)] =(x-l) k ^-\ v \T r (x,y). (3.8) 

Leaving aside the Tutte polynomial, several graph polynomials have been defined 
and intensively studied in the literature. The Tutte polynomial is a two-variable poly- 
nomial with one-variable specializations such as the chromatic polynomial or the flow 
polynomial. 

The chromatic polynomial is a graph polynomial Pr{k) (k G N*) which counts the 
number of distinct ways to color the graph T with k or fewer colors, where colorings are 
considered as distinct even if they differ only by permutation of colors. For a connected 
graph, this polynomial is related to the Tutte polynomial (13. 3p by the relation 

P r (k) = (-l)l^l- 1 A;T r (l - k, 0). (3.9) 

In order to define the flow polynomial, we need a finite Abelian group G. One can 
choose arbitrarily an orientation for each edge of the graph T, the result being inde- 
pendent of this choice. (The same type of situation appears when computing Feynman 
integrals, see the next section.) A G-flow on V is a mapping 

ip-.E^G (3.10) 

that satisfies current conservation at each vertex. A G-flow on T is said to be nowhere- 
zero if ip{e) 7^ for all e. Let Fp(G) be the number of nowhere-zero G-flows on T. One 
can prove that that this number depends only on the order k of the group G; it can 
thus be written Fr(k) — it is the restriction to non-negative integers of a polynomial 
in k, the flow polynomial. One can show that 

Fr(Jfe) = (-l) |i?h|y|+1 Tr(0, 1 - k). (3.11) 

4. QFT, Feynman integrals and their parametric representation 

In this section we introduce a few notions of QFT, the latter being a general frame- 
work lying at the very heart of fundamental physics. QFTs give a quantum description 
of particles and interactions, which is naturally compatible with Einstein's theory of 
special relativity. As already mentioned in the Introduction, QFTs led to the Standard 
Model of elementary particle physics, which is one of the best experimentally tested 
theories (for example in huge particle accelerators like the CERN's LEP). QFT's math- 
ematical formalism also successfully applies to other branches of theoretical physics, 
like condensed matter or statistical physics. 

The interested reader may consult any of the very good textbooks on QFTs, such as 

or \m. 



The simplest field theoretical model is the $ 4 scalar model; it consists of a single type 
of field $( scalar field, 

$ : M 4 ->■ K, (4.1) 

where the target space IK is either ffi. or C. 

Note that we use here 1R 4 , which corresponds to a Euclidean metric. If one needs to 
work with a Minkowskian metric, then M 4 needs to be replaced by R 1 ' 3 (the time being 
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singled out from the three spatial dimensions). The number 4=1 + 3 corresponds to 
the fact that we work on a four- dimensional space-time. 

For a field theoretical model to be defined, one needs to write down an action, which 
is, from a mathematical point of view, a functional of the fields of the model. For the 
real $ 4 model, the action reads 



S[$(x)]= I d A x \ V^$(x)^$(x) + ^m 2 $(x) 2 + V int [$(x)] 

L h 1 — 1 



(4.2) 



d 
dxfj. ' 



and 



where the real parameter m represents the mass. Furthermore, we have 

Vi n t[$(x)] stands for the interacting potential, namely ^^(x) 4 (A being a real parameter, 
called the coupling constant). 

The formula ( 14. 2 j) is written in configuration space (or direct space). By applying the 
Fourier transform, one can write it in momentum space (which is actually the space 
where the elementary particle physics computations are usually performed): 



r 4 

S[$]= [ d 4 p \Y^p^p^ + l m & + V 1T _ 

^ 4 u=l 2 



(4.3) 



Let us now further analyze this action. One has a quadratic part in $ and a non- 
quadratic one. The quadratic part corresponds to the propagation (the free theory), 
while the non-quadratic part corresponds to the interaction. 

The functional integration is introduced as the product of integrals at each space point 
x (multiplied by some irrelevant normalization factor), TD<f>{x) := J^Yl x j d<&{x). Nev- 
ertheless, such an infinite product of Lebesgue measures is mathematically ill-defined. 
For a well-defined QFT measure, the interested reader may consult, for example, the 
review article [62] . 

The definition of the partition function is 

Z := I Z>$(aOe~ 5[ * (a0] . 



(4.4) 

The physical information of a theory is encoded in n-point functions (or correlation 
functions) which are defined by 



x 2 , ■ 



x N 



£ty(x)$(xi)$(s 2 ) ■ ■ ■ <5>(x N )e~ m 



[$(x 1 )$(x 2 )...$(x N )}. 



(4.5) 



In general, one is unable to find exact expressions for these correlation functions. 
In such cases, the tool used in theoretical physics is the perturbative expansion, i.e., an 
expansion of the exponential above in powers of A. The coefficients of such an expansion 
are sums of multiple integrals (see ( I4.2p and ( 14. 5p ); the number of these integrals grows 
rapidly with increasing order in perturbation theory (i.e., power of the coupling constant 
A). 

To these multiple integrals, called Feynman integrals, one associates Feynman graphs, 
which are very useful in the organization of the expansion coefficients. For the $ 4 model 
exhibited here, the graphs have degree four at each vertex; moreover one has internal 
and external edges (see, for example, Figure |3]). 
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Figure 3. A $ 4 Feynman graph, with four internal edges (ei, e 2 , e 3 , e 4 ) 
and four external edges (fi, f 2 , /3, /k). It is a graph at the third order in 
perturbation theory (it has three vertices). 

Each of these graphs comes with its combinatorial weight (which is highly non-trivial 
because of the non- labeling of the edges). Furthermore, the Feynman integrals can 
then be manipulated using various techniques of analytic combinatorics (the Mellin 
transform or the saddle-point approximation). 

In momentum space, one orients each edge e (internal or external) and associates 
to it some momentum p e . The contribution from any internal edge comes from the 
quadratic part of the action (14. 3[) and reads 



This contribution is called the propagator. 

The contribution of some vertex is given by the coupling constant A times a 8- 
function of conservation of the incoming/outgoing momenta at the respective vertex. 
If the momenta of the external edges are fixed (they correspond to physical states), the 
momenta of the internal edges are free and one has to integrate over them, leading to 
the respective Feynman integral. 

For the example of Figure [3j all this reads 



We leave it as an exercise to the interested reader to find out the orientation of the 
edges which was chosen for such a Feynman integral to occur. 

Because of the presence of the three ^-functions in (14.71) . the number of remaining 
integrals is equal to two, which is actually the number of independent cycles of the 
Feynman graph (this being a general result in QFT). 

At the end of the story, for the Feynman integral (14.71) all this leads to a logarithmic 
divergence in the high energy (\p\ — > 00) regime of the internal momenta (the so- 
called ultraviolet regime). The appearance of such divergences is in fact a very frequent 
phenomenon in QFT; it is the renormalization process which, when possible, takes care 
of these infinities in a highly non-trivial way (see, for example, the book [61J). 

Let us now give some insight in this process of renormalization in QFT. If the con- 
sidered QFT model is renormalizable, the graphs which lead to divergences should 
correspond to terms present in the action. In order to illustrate what we mean by 




(4.6) 




• S(Pfi + Vh ~ Vex ~ Pe 2 )8(p ei + Pe 3 ~ Pe t ~ Pu)8(Pe 2 ~ Pe 3 + Pe 4 ~ P / 3 ) • (4.7) 
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this, let us go back to the example of the $ 4 model. The graphs which lead to the 
various divergences of the model need to have two or four external edges. Thus, these 
divergences can be "cured" by an appropriate renormalization of the parameters of the 
action H4.2[) (for example, the mass m and the coupling constant A). For instance, a 
graph with four external edges (each of these four external edges being associated to 
a field $) corresponds to the renormalization of the coupling constant A. This comes 
from the fact that the coupling constant A is indeed the parameter by which the $ 4 
term in the action (14.21) is multiplied. 

Let us take as an example the graph T with four external edges in Figure HI In 




Figure 4. A four-point Feynman graph. It has two vertices, localized, 
in configuration space, in x and y. 



configuration space, the Feynman integral to renormalize reads (up to irrelevant nor- 
malization factors) 

0(r) = J dxdy(C(x, y)) 2 C(x, y a )C(x, y 2 )C(y, y 3 )C(y, y A ), (4.8) 

where we have denoted the four external points that the external edges hook to by 
2/1)2/2)2/3,2/4- Note that the propagators C(x,y) are given by the inverse Fourier trans- 
forms of the momentum space propagators (14. 6ft — they are associated to the quadratic 
part of the configuration space action (14. 2p . i.e., they are integral representations of heat 
kernel functions (by a slight abuse of notation we have denoted propagators in both 
configuration and momentum space by C). 

We now subtract a term corresponding to the configuration space region 

x ~ y. (4.9) 

The integral (14. 8p is divergent and can be rewritten, using an appropriate Taylor ex- 
pansion, as 

dxdy(C(x, y)) 2 C(x, y\)C(x, y 2 ) 

C(x, y 3 )C{x, y A ) + / dt{y - x) V {C(x + t(y - x),y 3 )C{x + t(y - x),y A )) 



(4.10) 

One can prove that the second term in the sum above — the renormalized Feynman 
integral — is finite (see, for example, the book [HI] for details), while the first is the one 
containing the divergence — it will be subtracted. This term is called local counterterm: 

r r 0(r) = J dx{${x)Y J dy(C(x,y)) 2 . (4.11) 

The adjective "local" comes from the fact that this expression is associated, as shown 
above, to the configuration space region (14. 9p . The operator r represents the first term 
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in the Taylor expansion ( I4.10p . Note also that in (14. lip we have explicitly considered the 
multiplication with the four fields — this divergent term can thus be reabsorbed, 
as already mentioned above, in a redefinition of the coupling constant A. Since the 
model (14. 2p is translation- invariant, the second integral in (14. lip is independent of x, 
and thus the counterterm (14. lip is of the same form as the interaction term in the 
action. We have thus shown that the subtraction of a local counterterm makes the 
renormalized Feynman integral finite. 

The locality of the counterterms can be seen as sending all of the external edges to 
the same point (here x), again, in configuration space. The phenomenon described here 
can also be exhibited in momentum space (see, for example, the review paper [62J for 
the same example treated in momentum space). 

Nevertheless, subtraction of the divergences can be more involved. One can have 
graphs with sub divergences. In order to introduce the general (Bogoliubov) subtraction 
operator for a Feynman graph, we first need the following definition. 

Definition 4.1. A Zimmermann forest J 7 of a Feynman graph T is a set of subgraphs 
of T such that 

7 n i = or (7 C i or 7' C 7), for all 7, 7' G 7. (4.12) 

Remark 4.2. The Zimmermann forest defined above is a notion distinct from that of a 
forest in graph theory (see Definition 13.2( 7)). 

The (Bogoliubov) subtraction operator reads 

R:=J2U(- T ^ ( 4 - 13 ) 

where the sum is over all Zimmermann forests T of superficially divergent subgraphs 
7 (including the empty forest). The superficially divergent graphs are, in the case of 
the $ 4 model, the graphs with two or four external edges (this result comes from the 
so-called power counting theorem). 

Applying this subtraction operator to a Feynman integral 0(r) extracts the diver- 
gences. The remaining part, the renormalized Feynman integral 

«Mr) = (4-i4) 

is finite. 

An example of a Feynman graph with sub divergences is given in Figure where the 
Feynman graph of Figure H] is a subgraph. Renormalizing this subdivergence leads to 
the renormalization of the whole graph. This can be explicitly seen from the following 
facts. The Zimmermann forests of this graph are 

0,{ 7 }- (4.15) 

The renormalized Feynman integral thus reads 

+ (r) = R(<j>(T)) = 0(r) - (r 7 (0( 7 )))0(r/ 7 ), (4.16) 

where we have used the general property r 7 0(r) = <p(T/j)r(j)('y), 7 being a subdiver- 
gence of the Feynman graph T. Note that we wrote 7/7' for the cograph obtained by 
shrinking the subgraph 7' inside 7 — shrinking a subgraph means to replace its internal 
structure by a point. 
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Figure 5. A Feynman graph, with six internal edges and with a subdi- 
vergence given by the subgraph made of the internal edges 1 and 2. 

This graph-by-graph renormalization method is known under the name of the Bogo- 
liubov-Parasiuk-Hepp-Zimmermann (BPHZ) renormalization scheme, and it is one of 
the most frequently used in QFT. 

For a general Feynman graph, it is also possible to have overlapping divergences, 
i.e., divergent subgraphs with a non-empty intersection. In this case, one can use a 
classification theorem of Zimmermann forests in order to prove finiteness of renormal- 
ized Feynman amplitudes. In this case, one needs a scale decomposition of the graph's 
propagators. This means that each propagator C is written as a sum of distinct contri- 
butions, each such contribution corresponding to a distinct energy scale (see the book 
[pTj for a detailed presentation of these issues). 

Let us also mention that, once the renormalization techniques are performed, the 
renormalized (and hence finite) integral leads to physical quantities which are measured 
in elementary particle collider experiments with an extremely high accuracy. 

In the case presented here, the $ 4 model, the mass m prevents the integral to be 
divergent also in the infrared regime (that is, for \p\ — > 0; see, for example equation (14. 7p . 
which is convergent in this regime). 

For the sake of completeness, let us also mention that, when computing Feynman 
integrals, mathematical physicists are interested in graphs with no bridges (called 1- 
particle irreducible (1PI) graphs). This comes from the fact that the momentum of the 
bridge cannot flow towards its high energy regime, because it is directly related to the 
external (and thus always finite) momenta. These graphs are thus of no interest for a 
renormalization analysis. 

Let us now proceed further and introduce the parametric representation of a Feynman 
integral. The main idea is to write each of the internal propagators (14. 6 p of the integral 
as an integral over some parameter a: 



Jo 

Inserting these formulas in the general expression of a Feynman amplitude allows to 
integrate out (through Gaussian integrations) the internal momenta p e . 




e = 



1,2,. 



(4.17) 
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The Feynman integral then becomes 

roo -VHp cxt ,a)/[/ r (a) W 

m = L w ffi"" °' da ' ) - (4 ' 18) 

where t/r(«) is a polynomial in the set of parameters a and Vr(p ex t, «) is a polynomial 
in the set of external momenta p ex t and the set of parameters a. It can be shown that 
these polynomials depend only on the underlying graph. More specifically, one has 

f/ r(a) = ^n^ ( 419 ) 
r e&r 

and 2 

Vr(p ex t,a) = J2H^( Yl Pi) > ( 42 °) 
r 2 tph. Se£(75) ' 

where the first sum is over all trees T of the graph T, while the second sum is over 
two-trees %, which, as already stated in the previous section (see Definition 13.2( 9)) 
separates the graph in two connected components and one can identify the external 
edges which connect to one or the other of these connected components. The symbol 
E(J~2) denotes one of the connected components thus obtained. 

Remark 4.3. By momentum conservation, the total momentum of one of these connected 
components (for example E(7^)) is equal to the total momentum of the other connected 
component. 

For the example of Figure [31 one has 

U-p(ot) = a^a^ + ol-iOl^ + a^c^ + «i«3 + «i«4, 
Vr(Pext,a) = (p/i + Pf 2 ) 2 aia 2 (a 3 + a 4 ) + p 2 f4 aia 3 a 4 + p 2 f3 a 2 a 3 a>4. 

5. Relation between the multivariate Tutte polynomial and the 
polynomials of the parametric representation 

Let us now exhibit the relation between the notions presented in the two previous 
sections, the multivariate Tutte polynomial on the one hand, and the polynomials of 
the parametric representation of Feynman integrals on the other hand. 

We start with a deletion/contraction relation for the polynomials Ur(a). 

Theorem 5.1. For any semi-regular edge e, we have 

U T (a) = a e U T _ e (a) + U T/e (a). (5.1) 
Moreover, the terminal form evaluation is 

U r (a) = Y[a e , (5.2) 

e 

for G consisting only of self-loops attached to isolated vertices. 

In order to prove this result, we shall use the Grassmann representation (see Appen- 
dix [SJ of the polynomial £/r. 

A key ingredient for this representation is the incidence matrix (e ejV ) of the graph. 
In this matrix, the rows are indexed by the edges e of the graph, e = 1, 2, ... E, and 
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the columns are indexed by the vertices v of the graph, v = 1, 2, . . . , V. By definition, 
the entry e e>v is equal to 1 if the edge e is outgoing from the vertex v, — 1 if the edge e 
is incoming at the vertex v, and it is if the edge e is not incident to the vertex v. Let 
us point out that this matrix encodes the complete information of a Feynman graph 
without self-loops. In order to encode the information from self-loops also (if they are 
present), one needs to define an additional matrix (i]e,v), with non- vanishing entries r] ejV 
if e is a self-loop incident to the vertex v and r\ e<u = otherwise (note that e e> „ = if e 
is a self-loop incident to the vertex v, see [44J for details). 

The proof of Theorem 15. II is then based on the following auxiliary result. 

Lemma 5.2. The polynomial Ur(a) is given by 

detQr, (5.3) 

where the graph matrix Qr is the square matrix 

*=(■(-;,) '"o"')- < 5 - 4) 

{Here, *( • ) denotes the transpose of a matrix) . 

Proof. One can prove this by performing the Gaussian integrations obtained from in- 
serting the propagator form (I4.17P in the general form of a Feynman integral (see again 
for details). □ 



Proof of Theorem \5.1[ As a direct consequence of Lemmas 15.21 and IA.21 we have 

U G {a e )= I \\d Xv du v dXedu e e- aeX ^e- x ^ xMe ^. (5.5) 

v,e 

Note that we have split each of the two sets of Grassmann variables (x and u) required 
for the representation of such a determinant (see Lemma IA.2j) in two subsets, one of 
cardinality E (indexed by e) and one of cardinality V (indexed by v); this corresponds 
to the block form of the complete matrix Qr- 

Let now e be a semi- regular edge connecting vertices V\ and v%. We have (see the 
definition (lA.4j) for a single Grassmann variable) 

e" aeXe " e = 1 + a e u eXe . (5.6) 

We denote the two terms of this sum by det Qr, e ,i and by detQr,e,2, respectively. Let 
us first closely investigate the term det Qr,e,i- By direct inspection of the various terms 
coming from similar expansions of the other Grassmann exponentials in (15.51) . one can 
prove (using the Grassmann integration rules (IA.5I) ) that (see again [44] for details) 



det 



Qr,e,i = j \\ dxe'du e >dxvduj v (xv 1 - Xv 2 ){^v 1 - 



Let us now perform the following (triangular) change of variables (with unit Jacobian) 

Xvi Xv\ Xv2 ) 

Xv = Xv, for all v ^ v u 
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and the same for the set of u Grassmann variables. This leads to the identification (see 
again [H] for details) 

det Qr, e ,i = det Q r / e - (5.8) 
The second term of the expression (15.61) similarly leads to 

det Qr, e ,2 = « e det Qr- e - (5.9) 

We refer the interested reader again to [H] for obtaining the terminal form evaluation. 

□ 

Remark 5.3. The Grassmann calculus method presented here for the proof of Theo- 
rem 15.11 is not the simplest method to obtain this result (which can be obtained for 
example using the explicit form ( I4.19P of the polynomial Ur{ct)). Nevertheless, the 
Grassmann calculus method is particularly useful to get, for example, formula ( I4.19P 
(or formula (14.201) for the polynomial Vr(p ex t, ot)) directly from the respective Feynman 
integral. Moreover, such a method can be used to obtain the parametric representation 
of more complicated QFT models, like for example the non-commutative scalar <3>* 4 
model (see Subsection 7.2). 



The situation is similar for the second polynomial V(p ext , ex) (see [2] for details). Let 
us now show how the polynomial Ur{ot) can be obtained as a limit of the multivariate 
Tutte polynomial Z G (q,(3). One needs to consider 

q- k ^Z T {q^). (5.10) 

Taking the limit q — > one obtains a sum over maximally spanning subgraphs A, that 
is, subgraphs with k(A) = k(T): 

S r (P)=\imq~ k ^Z r (q,P)= V Y[/3 e (5.11) 

A maximally spanning E e£A 

If one now retains only the terms of lowest degree of homogeneity in /3, one obtains 
a sum over maximally spanning graphs with lowest number of edges, i.e., maximally 
spanning acyclic graphs (or spanning forests, see Definition 13.2( 7)) of T. We denote 
this lowest number of edges by p (note that p = \ V(T) \ — k(T)). We have 

T maximally spanning forest of G edT 

Finally, if one divides Fr((3) by Y\ e&E (3 e and makes the change of variables a e = 
the polynomial U^{a) is obtained. 

All this is summarized by the formula 



U r (a) 



nf3 e ) lim -=r lim q- k{r) Z r (q, q'(3) 



(5.13) 



Note that the variable q' was introduced to take the limit retaining only the lowest 
degree of homogeneity in (3 from Sr(q'P) (see above). 

Let us also mention here the universality of the Tutte polynomial. This, also known 
as the "recipe theorem", means that any Tutte-Grothendieck invariant must be an 
evaluation of the Tutte polynomial, with the necessary substitutions given by the recipe 
(see, for example, the survey papers [281 |29]). The relation between the Symanzik 
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polynomial Ur and the multivariate Tutte polynomial can thus be deduced from 
this result and Theorem 15.11 

6. Combinatorial Connes-Kreimer Hopf algebra of Feynman graphs 

Consider now the unital associative algebra H freely generated by 1PI Feynman 
graphs of the $ 4 model (including the empty set, which we denote by 1). The product 
m is bilinear, commutative, and given by the operation of disjoint union. 

Let us give a formal definition of the Feynman rules, already introduced in the pre- 
vious section. 

Definition 6.1. The Feynman rules are a homomorphism from H to some target 
space A. 

For the sake of completeness, let us also remark that this target algebra is naturally 
equipped with some Rota-Baxter algebraic structure. Recently, the algebraic implica- 
tions one has when relaxing this particular condition were studied in [27]. 

In a generic QFT model, a special role in the process of renormalization is played by 
the primitively divergent graphs. 

Definition 6.2. A primitively divergent graph of a QFT model is a graph whose Feyn- 
man integral is divergent but which does not contain any subgraph for which the Feyn- 
man integral is also divergent. 

In the particular case of the commutative $ 4 model, this class of graphs is formed 
by the graphs with two or four external edges that do not contain any graph of two or 
four external edges as subgraph. Let us also remark that the parametric representation 
presented in Section @] is one of the elegant ways to prove this type of result. 

Definition 6.3. The projection T is a map from A to A which satisfies the following 
property: For all T G "H, T primitively divergent, 



This means that if 4>{T) is divergent then its overall divergence is completely included 
inTo0(r). 

From now on we denote the set of two- and four-external edges subgraphs of a generic 
graph T by T. This is also referred to as the set of superficially divergent subgraphs of 

r. 

Lemma 6.4 ( (TTJ Lemma 3.2]). Let T G "H. Provided 

(1) for all 7 G T_ and 7' G 7, the graph 7/7' is superficially divergent, 

(2) for all 71 G H and 72 G "H such that 71 and 72 are superficially divergent, there 
exists gluing data G such that (7! o G 72) is superficially divergent, 

the following coproduct is coassociative: 



(\& A - T) o 4>(T) < 00. 



(6.1) 



Ar = r®i + i®r + at, 



(6.2a) 
(6.2b) 
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By gluing data we understand a bijection between the external edges of the graph to 
be inserted and the edges of the propagator (or vertex) where the insertion is done. 

Let us now illustrate how this result fits the commutative $ 4 model. The first condi- 
tion of the lemma is trivial, since when one shrinks a subgraph, the number of external 
edges is conserved. Moreover, since 7' has two or four external edges, the cograph is 
also a $ 4 graph. Let us now check condition [5] of Lemma 16.41 We consider two graphs 
71 and 72 with two or four external edges. We consider 70 = 71 o G 72 for any gluing 
data G. Let Vj, U and Ei be the numbers of vertices, internal, and external edges of 
7t, i G {0, 1, 2}, respectively. For all i G {0, 1, 2}, we have 

4V- = 2Ii + E h (6.3a) 
\ A 1 ' Vj - (6.3b) 



V 1 + V 2 -l, 

h + h + l, 
h + h, 



if E 2 


II 


if E 2 


= 4, 


if E 2 ■- 


= 2, 


if E 2 


= 4, 



(6.3c) 



which proves that E = Ei. Then, as soon as 71 is primitively divergent, so is 70. 
Concerning condition [TJ note that 7" = 7/7' if and only if there exists G such that 
j = y 0(5 7', which allows to prove that condition [TJ also holds and that the coproduct 
(16. 2 p is coassociative. 

Let the coproduct A : H — > H <8> % be defined by 

Ar = r®i + i®r + ^7®r/ 7 , vre-H. (6.4) 

7 er 

Furthermore, let us define the counit e : H — > K by 

e(l) = 1, e(r) = 0, for all T ^ 1. (6.5) 

This means that, for any non-trivial element of H, the counit returns a trivial answer 
and, equivalently, only for the trivial element of H (the empty graph 1), the result 
returned by the counit is non-trivial. Finally the antipode is given recursively by 

S(l) = I, 

T^-T-J2S(l)Th. (6.6) 

Then the following result holds. 
Theorem 6.5 ([2TJ Theorem 1]). The quadruple (H, A, e, S) is a Hopf algebra. 
Let us also notice that Tl is graded by the loop number. 

Let now f,g G Hom("H,^4), where A is the range algebra of the projection T (see 
above). The convolution product * in Hom("H,^4) is defined by 

f*g = m A o(f®g)oA n . (6.7) 
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Let be the Feynman rules and 0_ G Hom('H, A) the twisted antipode, defined recur- 
sively by 

0_(r) = -r(0(r) + ^0_( T )0(r/ T )), foraiireH. (6.8) 

Let us mention here that this twisted antipode implements the BPHZ subtraction (see 
Section 4). The renormalized Feynman amplitude (f) + of a graph r e % is given by 

+ (r) = 0_*0(r). (6.9) 

Moreover, Bogoliubov's subtraction operator is given by (see, for example, |37] for 
details) 

R(T) = 0(r) + 0-(7)0(r/ 7 )- (6.10) 

Let us now illustrate this with the example of the divergent Feynman graph of Figure [5j 
In that case, the expression (I6.10P above becomes 

J R(0(r)) = 0(r) + 0_( 7 )0(r/ 7 ), (e.ii) 

where the only sub divergent graph 7 is given by the internal edges 1 and 2 (see again 
Figure |5]). The above definition of the twisted antipode leads to 

<M 7 ) = -T(0( 7 )) (6.12) 

(since the subgraph 7 is primitive — the recursion in the definition of the twisted 
antipode stops). Inserting the expression (16.121) for the twisted antipode in (16. lip 
finally gives 

R(<f>(T)) = 0(r) - T(0( 7 ))0(r/ 7 ), (6.13) 

which is, as announced, the same expression as (14. 16j) . the expression computed using 
non-algebraic formulas. Note that the explicit form of the projection T is given by the 
choice of the renormalization scheme. In Section 4, working in configuration space, we 
have used explicit Taylor expansion operators r. 

For the sake of completeness, let us recall that a combinatorial Hopf algebraic struc- 
ture like the one described in this section can be also defined for more involved quantum 
field theories, such as gauge theories [4"B"1 [55]. 

Moreover, let us also state that a slightly different combinatorial Hopf algebra — the 
core Hopf algebra — was defined in [S] (and independently in [13] for vacuum graphs). 
Its coproduct does not sum over the class of superficially divergent subgraphs but over 
all subgraphs of the respective graph: 

Ar = r<g>i + i<g>r + ^ 7 <g>r/ 7 , vr en. (6.14) 

The definitions of the coproduct and of the counit follow in a straightforward manner. 
One can then verify that this structure is also a Hopf algebra structure. The cohomology 
of this Hopf algebra was also investigated in detail [4"8l 139]. 

This core algebra contains the renormalization Hopf algebra (presented in this sec- 
tion) as a quotient algebra. Moreover, it can be seen as the Connes-Kreimer Hopf 
algebra of a QFT model formulated in infinite dimension (see again [4"8| |4"9] for more 
details). 
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Let us end this section by mentioning that the study of the cohomology of the Connes- 
Kreimer Hopf algebras allows one to express the combinatorial Dyson-Schwinger equa- 
tions as power series in some appropriate insertion operators of graphs (which, from a 
mathematical point of view, are Hochschild one-cocyles). The interested reader may 
consult the articles H6j [50] or the thesis (75]. For various results on these combina- 
torial Dyson-Schwinger equations in the Hopf algebra of decorated rooted trees, one 
may also consult [301 EI] • 

7. Ribbon graphs and scalar QFT on the non-commutative Moyal space 

In this section we present the generalization of the results of the previous sec- 
tions to the case of ribbon graphs (where the Tutte polynomial generalizes to the 
Bollobas-Riordan polynomial) and the case of $ 4 scalar QFT on the 4-dimensional 
non- commutative Moyal space. 

7.1. Ribbon graphs — the Bollobas-Riordan polynomial. One has the following 
definition, generalizing in a natural manner the definition of a graph. 

Definition 7.1. A ribbon graph F is an orientable surface with boundary represented as 
the union of closed disks, also called vertices, and ribbons also called edges, such that: 
the disks and the ribbons intersect in disjoint line segments, each such line segment lies 
on the boundary of precisely one disk and one ribbon and, finally, every ribbon contains 
two such line segments. 

Let us also mention that ribbon graphs can also be defined as graphs equipped with, 
for each vertex, a cyclic ordering of the edges incident to that vertex, or as graphs 
embedded in surfaces (this last definition was actually the object on which B. Bollobas 
and O. Riordan defined their generalization of the Tutte polynomial in [5| [TP]). 

The definition above can be extended such that a distinct type of edge — the external 
edge — is permitted. This type of edge contains only one line segment given by the 
intersection between a ribbon and a disk (see Definition 17. ip . 

An example of such a graph is given in Figure [6j 




Figure 6. An example of a ribbon graph with one vertex, one internal 
edge, and two external edges. 

Definition 7.2. A face of a ribbon graph is a connected component of its boundary as 
a surface. 

For example, the graph of Figure [6] has two faces. 

If we glue disks along the faces, we obtain a closed Riemann surface whose genus is 
also called the genus of the graph. 

Definition 7.3. The ribbon graph is called planar if this Riemann surface has genus 
zero. 

For example, the graph of Figure [6] is planar. 
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Definition 7.4. A planar graph is called regular if the number of faces broken by 
external edges is equal to 1. 

The graph of Figure [6] is planar irregular. 

As already stated in the Introduction, ribbon graphs are also known as maps (see 
for example [18] and references therein). Planar regular maps are also known as outer 
maps. 

Definition 7.5. The Bollobas-Riordan polynomial of a ribbon graph G is defined by 
R T (x,y,z) := J2(x-l) rir) - r{H V {H) z HH) - FiH)+niH) - (7-1) 

HcE 

In the above definition, the symbol F(H) stands for the number of components of 
the boundary of the respective subgraph H (the number of faces). The supplementary 
variable z is required to keep track of the topological information (the graph genus or 
the number of faces). 

As the Tutte polynomial, the Bollobas-Riordan polynomial also obeys a deletion/con- 
traction relation. 

Theorem 7.6. Let V be a ribbon graph. Then 

R r = Rr/e + Rr-e (7.2) 

for any regular edge e of G, and 

R r = xR r/e (7.3) 

for every bridge of T. 

There are further analogies to the Tutte polynomial, in the sense that, defining 
the Bollobas-Riordan polynomial on terminal forms transforms property (17. 2p into a 
definition. On these terminal forms (that is, graphs with one vertex), the polynomial 
is defined by 

R r (y,z):=^y lEm ^ 9iH \ (7-4) 
Her 

since, in this case, one has k(H) — F(H) + n(H) = 2g(H). 

A multivariate version of the Bollobas-Riordan polynomial exists as well in the lit- 
erature, namely 

Z r (x,f3,z)=J2* k(H) (ll( 3 ) zFiH) - ( 7 - 5 ) 

HcE \eeH / 

This version also satisfies a deletion/contraction relation. 

Finally, let us mention that a signed version of the Bollobas-Riordan was defined in 
[20J; this is a three- variable polynomial defined on signed ribbon graphs, that is, ribbon 
graphs for which an element of the set {+, — } is assigned to each edge. A partial duality 
with respect to a spanning subgraph was also defined [!§]; this allows one to prove that 
the Kauffman bracket of a virtual link diagram is equal to the signed Bollobas-Riordan 
polynomial of some ribbon graph constructed from a state of the respective virtual link 
diagram (the interested reader is referred to [19] for details on this topic). Moreover, 
the properties of the multivariate version of this signed Bollobas-Riordan polynomial 
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were analyzed in [73] (namely, there its invariance under the partial duality of [19] 
was proved) . Finally, 4- variable generalizations of the Bollobas-Riordan polynomial for 
ribbon graphs were defined in [51] and [33]. 

7.2. Scalar $* 4 QFT on the non-commutative Moyal space; parametric rep- 
resentation. Let us start this subsection with some explanations on the interest of 
mathematical physicists in QFT on non- commutative spaces. 

Non-commutative QFT can be seen as a potential candidate for new physics be- 
yond the celebrated Standard Model of elementary particle physics. Moreover, non- 
commutativity of space-time is believed to be a possible framework for the quantification 
of gravity. Thus, it was argued in [25] that the space-time becomes non-commutative at 
the Planck scale (the physical energy scale where gravitational effects have the same or- 
der of magnitude as quantum effects; thus, the energy scale where gravity and quantum 
physics cannot "ignore themselves" anymore). Moreover, non-commutative QFT has 
been shown to be related to loop quantum gravity and to group field theory [3j 134"! 133]. 
some of the major candidates for a quantified theory of gravity (more details on group 
field theory are also given in the last section of this survey). Non- commutative mod- 
els have also been shown to be effective theories of certain models of string theories 
and of matrix theories [221 E3]- Finally, changing point of view, it was argued that 
non-commutative models can be better adapted for describing the physics of effective 
non-local interactions, as it is the case in the fractional quantum Hall effect [331 EO] • For 
more details on non-commutative QFT, the interested reader may consult the review 
paper [63] . 

From a mathematical point of view, the Moyal space, one of the simplest cases of 
such a non- commutative setting, can be defined in the following way. 

Definition 7.7. The Moyal algebra is the linear space of smooth and rapidly decreasing 
functions iS(M' D ) equipped with the Moyal product 

(f*g)(x) = J -^d D yf{x+ l -Q-k)g{x + y)e lk -y^xeM D , (7.6) 

where G is a D-dimensional skew-symmetric matrix. 
Remark 7.8. 

(1) The Moyal product is non-local, non-commutative, but it is still associative. 

(2) Taking the matrix to be vanishing leads to the usual commutative multipli- 
cation of functions. 

The Moyal algebra defined above can be extended by duality, considering the product 
of a tempered distribution with a Schwartz-class function. The identity, the polynomi- 
als, the 5 function, and its derivatives can then be shown to belong to this extended 
Moyal algebra (see [35] for details). One can then prove that 

[xp, Xy]* := x^ -k x u — x u -k x^ = iO M ", for all fi, v — 1, 2, . . . , D. (7.7) 

Let us note that this relation can also be used as a definition of the Moyal space. 
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In this survey we take D = 4 (the usual 4-dimensional space-time). The non- 
commutativity matrix 9 is taken to be 

e - f 02 M e - f° ~ e 
~ V o e 2 J ' ° 2 " [e o 

Non-commutative quantum field theories are then implemented by replacing the usual 
commutative multiplication of fields by the non-commutative ^-product. Thus, the 
action of the non- commutative $* 4 model reads 

4 



(7.8) 



where the non-commutative interaction potential is given by 

km = ^* 4 (*)- (7-9) 

As in the commutative case, one can apply the Fourier transform in order to transfer 
this action to one on momentum space. A first consequence of the use of the Moyal 
product is that its use does not change the propagation part of the action. This is a 
direct consequence of the following mathematical property: 

y"d 4 £($*^)(£) = J d 4 x$(x)^(x), for all $ and ^. (7.10) 

Thus, the propagator of the theory has the same form as in the commutative case; this 
is particularly useful for implementing the parametric representation in a similar way. 

Nevertheless, the use of the Moyal product at the level of the interaction terms 
changes completely the situation: the non-commutative interaction is not equivalent to 
the commutative one (as this is also the case for the propagation). The interaction part 
no longer preserves the invariance under permutation of external fields. This invariance 
is restricted to cyclic permutations only. Furthermore, there exists a basis — the matrix 
base — of the Moyal algebra where the Moyal product takes the form of an ordinary 
matrix product. For these reasons the associated Feynman graphs are ribbon graphs, 
just like the ones described in the previous subsection. 

For more details on the implementation of QFTs on the non-commutative Moyal 
space, the interested reader is referred to [H3] and references therein. 



Let us now move on to the implementation of the non- commutative parametric rep- 
resentation. This can be achieved using the same trick (I4.17P as in the commutative 
case. The difference comes however from the fact that the interaction part is non- 
trivially dependent on the non-commutativity parameter 9, and this dependence will 
consequently appear in the form of the two polynomials which give the parametric rep- 
resentation of the respective Feynman integral (we denote these polynomials by U*(a, 9) 
and V*(p cxt ,a,9)). 

Using a Grassmannian development of the Pfaflians obtained through these Gaussian 
integrations, one can prove (see again [44J) that the Feynman integral of a generic graph 
T reads 

poo -V*(pext,a,6)/U*(a) & 
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In order to present the exact form of the polynomials U£(a,8) and V^(p ext , a, 9),we 
need the following definition (see again [H] for details). 

Definition 7.9. 

(1) A *-tree of a connected ribbon graph is a sub-ribbon graph with one boundary 
component. 

(2) A two *-tree is a sub-ribbon graph with two boundary components. 

Note that these definitions are non-trivial generalizations of the notions of trees and 
two-trees, respectively (see Definition 13.2( 8) and (9)). Let us also mention that *-trees 
are also known as quasi-trees in the knot theory literature and as unicellular maps in 
the combinatorics literature (see for example [73] and [17] and references therein). 

We have 




(7.12) 



where b = F — 1 + 2g. 

The case of the polynomial V*(p ext ,a,6) is more involved, because the non-com- 
mutative setting presented here allows for it to be complex. We have the following 
theorems. 

Theorem 7.10 ([HJ Theorem 5.2]). The real part of the polynomial Vp* {p e xt, ot, 9) reads 
X£(p ext ,a,9)= (£) 6+1 £ II 2 y^) 2 , (7-13) 

7^* two-*-tree e^T^ 

where p-r* is the sum of the external momenta entering one of the two faces of the 
two-k-tree T£ ■ 

Note that we say that an external momentum enters a face if the respective external 
edge is hooked to one of the vertices of the face (we also consider in the sum above, 
for each such external momentum, the sign + if the momentum is ingoing, respectively 
— if it is outgoing). As in the commutative case, the choice of the face in the above 
theorem is irrelevant (by momentum conservation). 

Theorem 7.11 ((HI Theorem 5.3]). The imaginary part of the polynomial Vp (p e xt> ot, 8) 
reads 

3$(p«t,M) = (l) E U 2 T^ ( 7 - 14 ) 

where ip(p) is the phase obtained by following the momenta entering the face of the 
•k-tree T*. 

Let us give some explanations on the definition of ip{p). We denote the (unique) face 
of the *-tree by /. One can then define an incidence matrix (e/,i) between this face and 
the external edges of the graph (in the same way one defines the usual incidence matrix 
between vertices and edges). Following the external momenta entering the face, one can 
define an order relation on these momenta. One then has ip(p) = ^2i<ji^f,iPi) ^ {^f,iPj)i 
where Pi Apj = \p^Q^ v 'P l j- This expression generalizes in a straightforward manner the 
expression which gives, as a function of the incoming/outgoing momenta, the phase 
corresponding to a Moyal vertex. 
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7.3. Relation between the multivariate Bollobas— Riordan polynomial and the 
polynomials of the parametric representation. It can be shown that also the 
polynomials U£(a,9)r and V^(p ext ,a,9) satisfy a deletion/contraction relation. More 
precisely, we have the following result. 

Theorem 7.12 ([HI Theorem 5.1]). For any semi-regular edge e, we have 

U*(a, 6) = a e ll*_ e (a, 9) + U£ /e (a, 9). (7.15) 

One can moreover prove the following relation between the multivariate Bollobas- 
Riordan polynomial and the polynomial U£(a,9) of the parametric representation: 

U*(a, 9) = {9/2)W\~\^) ( \{ a.) x lim w~ x Z v (l, w) . (7.16) 

e€E 

Note that the commutative limit (6 — > 0) gives the usual $ 4 model on the commutative 
M 4 . This can be seen as well in the parametric representation, where U£(a, 9) is given by 
U r (a) plus some noncommutative corrections which factor 9 2 (corrections thus vanishing 
in the commutative limit). This directly applies to equation (I7.16P (where one obviously 
has to take the limits in the appropriate order on the right-hand side), and it leads to 

U T {a) = Mmt£(M) = lim ((9/2)^^)^^ x Hm«,-%(l, . 

(7.17) 

The second polynomial, V^(a,p ex t,9), can also be related to the multivariate Bol- 
lobas-Riordan polynomial in a similar way (see again jH] for details). 

For the sake of completeness, let us also mention that the parametric representa- 
tion of a non- commutative renormalizable scalar model, different from (17. 8p (the latter 
being actually non-renormalizable) — the Grosse-Wulkenhaar model [3B] — was fully 
implemented. Its polynomials do not obey a deletion/contraction relation, but they 
obey some other types of recursive relations [35]. Algebraic topological properties of 
the parametric representation of the Grosse-Wulkenhaar model (as well as of the para- 
metric representation of the commutative model (14. 2p ) were studied in [2]. 

8. COMBINATORIAL CONNES-KREIMER HOPF ALGEBRAS FOR RENORMALIZABLE 
$* 4 QFT ON THE NON-COMMUTATIVE MOYAL SPACE 

The definition of the Hopf algebra of non- commutative Feynman graphs which drives 
the combinatorics of renormalization is formally the same as in the commutative case. 
Nevertheless, in the definition of the coproduct one has to sum over the superficially 
divergent subgraphs of the respective graph, and this class of superficially divergent 
graphs takes now explicitly into account the topology of the graph. Thus, using var- 
ious QFT techniques, it was proved in [3B] that the superficially divergent graphs of 
the renormalizable Grosse-Wulkenhaar non-commutative model are the planar regular 
ribbon graphs with two and four external edges. 

We denote by W the unital, associative algebra freely generated by 1PI ribbon 
Feynman graphs of the renormalizable Grosse-Wulkenhaar $* 4 model (including the 
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empty set, which we denote by 1%*). The product m is bilinear, commutative and given 
by the operation of disjoint union. Let the coproduct A : W — > H* <S> H* be defined by 

Ar = l n * + In* ®r + ^ 7 ®r/ 7 , foralirGH*, (8.1) 

where, as earlier, we write £ for the set of two and four external edges planar regular 
sub-ribbon graphs of T. Furthermore, we define the counit e : H* — > K by 

e(l w *) = 1, e(r) = 0, Vr ± l w . (8.2) 

Finally, the ant ip ode is given recursively by 

S : H* -^H* (8.3) 

7 er 

Then the following result holds true. 

Theorem 8.1 f [7T| Theorem 4.1]). The quadruple (7i*, A, e, S) is a Hopf algebra. 

The main difficulty in the proof of this theorem (in comparison with the original 
Connes-Kreimer theorem) comes from the necessity of defining gluing data which re- 
spect to the ribbon graph cyclic ordering. Figure [7] shows an example of an insertion 
of a graph with four external edges with gluing data which respects the cyclic ordering 
(the dotted lines showing how the four external legs of the subgraph to be inserted and 
the four edges of the chosen vertex are put together — the insertion place). 




Figure 7. Insertion of a regular ribbon graph with four external edges. 
The cyclic ordering is respected by the chosen gluing data. 



The renormalized Feynman amplitude is defined in analogy to the commutative case 
(see equation ( 16. 91) ). 

Let us end this section by recalling that in [TO] the role played by Hochschild coho- 
mology of these types of Hopf algebras in non- commutative QFT was studied in detail; 
non-trivial examples of graphs with one or two independent cycles have been explicitly 
worked out. 
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Figure 8. A possible vertex of a three-dimensional quantum gravity 
tensor model. 

12 3 4 




10 2 f) 9 9 8 3 7 

Figure 9. A possible vertex of a four-dimensional quantum gravity ten- 
sor model. 

9. Perspectives — combinatorics of quantum gravity tensor models 

One of the main perspectives of the combinatorial approaches presented here is their 
extension for the study of the properties of quantum gravity models. The group field 
theory formalism of quantum gravity (for general reviews, see [32j [561 E3 EE]) is the 
most adapted for such a study, since it is formulated as a QFT. These models were 
developed as a generalization of two-dimensional matrix models (which naturally make 
use of ribbon graphs just as non-commutative quantum field theories do) to the three- 
and four-dimensional cases. 

The natural candidates for generalizations of matrix models in higher dimensions 
(> 2) are tensor models. In the combinatorially simplest case, the elementary cells that, 
by gluing together form the space itself, are the D-simplices (D being the dimension 
of space). Since a .D-simplex has D + 1 facets on its boundary, the backbone of group 
field theoretical models in dimension D should be some abstract interaction on 

rank D tensor fields <fi (see Figures [8] and [91 which represent, for D = 3 and D = 4, the 
vertex of these Feynman quantum gravity tensor graphs). 

In [37J ES] , some propositions for generalizations of the Bollobas-Riordan polynomial 
to the level of these rank three tensor models have been established, and it was proven 
that they respect a deletion/contraction property. Due to the increased topological 
complexity of the graphs, a different way of contracting edges had to be used. 

Furthermore, a polynomial for triangulations (no graphs being taken into consider- 
ations) has also been recently proposed in [52]. It is thus interesting to understand if 
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some connection between these propositions exists and if any of them could be related 
to some parametric representation of these quantum gravity models, as shown here to 
be the case for both commutative and non-commutative QFTs. 

Let us also stress the fact that different topological and analytic insights in this type 
of formulation of quantum gravity models have been investigated in the recent literature 
(see 01 El El El S21 E31 E2] and references therein) . 

It would be interesting to investigate whether or not these models can be related to 
the operation of cabling in knot theory and thereby to the colored Jones polynomials. 
Within the quantum gravity framework described in this section, (edge-) colored tensor 
models have been recently proposed and investigated (see [38] and references therein). 
On the other hand, the relation of parallel cabling to combinatorics in general and to 
the Bollobas-Riordan polynomial of a ribbon graph in particular was analyzed in |39j. 

For the sake of completeness, let us end this survey by mentioning other types of 
mathematical perspectives, such as number theoretical conjectures related to Feynman 
integrals, definition of knot invariants, etc. The interested reader may for example 
consult [13] and [T2~| [HI HSl EES] (and references therein) for various other developments 
related to the combinatorics of graph polynomials and QFT. 
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Appendix A. Grassmann variables, determinants and Pfaffians 

Grassmann variables Xi>X2> • • • , Xn are defined through their anticommutation rela- 
tions 

XiXj = -XjXh for all i and j. (A.4) 

A direct consequence of this relation is that any function of these variables is a poly- 
nomial with highest degree one in each variable. The rules of Grassmann integrations 
are then 

I dx=o and M= l (A - 5) 

The determinant of any n- dimensional square matrix can be expressed as a Grass- 
mann Gaussian integral over 2n independent Grassmann variables which one can de- 
note by V>2> • • • j ^ni ipii "02, • • • > VVi- (The reader should be warned that the bars have 
nothing to do with complex conjugation.) More precisely, we have 

detM = J JJ^d^e -E «* Mi ^. (A.6) 

The Pfaffian Pi(A) of a skew-symmetric matrix A is defined by 

detA= [Pf(A)] 2 , (A.7) 
where the sign of the term Ai^A^a • • • is fixed to be +1. 
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Proposition A.l ([1, Proposition 3]). The Pfaffian of an n- dimensional skew-sym- 
metric matrix A is given by 

Pf(A) = jY[dxie~^ XlA ^ = J dxidx2---dXne~^^ XlAl]Xj . (A.8) 

Lemma A. 2 ( |44| Lemma 2.1]). The determinant of a matrix D + A, where D is 
diagonal and A skew- symmetric, can be written in the form 

det(D + A) = f]J d Xl du l e- E * » B « W *-Si<i wAw+E^ Ui Ai jUj ( Ag ) 
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